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The genetic origins of Japanese populations have been controversial. Upper PaleoUthic Japanese, i.e. Jomon, 
developed independently in Japanese islands for more than 10,000 years until the isolation was ended with 
the influxes of continental immigrants about 2,000 years ago. However, the knowledge of origin of Jomon 
and its contribution to the genetic pool of contemporary Japanese is still limited, albeit the extensive studies 
using mtDNA and Y chromosomes. In this report, we aimed to infer the origin of Jomon and to estimate its 
contribution to Japanese by fitting an admixture model with missing data from Jomon to a genome-wide 
data from 94 worldwide populations. Our results showed that the genetic contributions of Jomon, the 
Paleolithic contingent in Japanese, are 54.3~62.3% in Ryukyuans and 23.1~39.5% in mainland Japanese, 
respectively. Utilizing inferred allele frequencies of the Jomon population, we further showed the Paleolithic 
contingent in Japanese had a Northeast Asia origin. 

The genetic origins of modern Japanese have been debated and three models are noteworthy'. In the 'con- 
tinuity' model, modern Japanese are considered as direct decedents of Jomon, the inhabitants of Japan in 
Paleolithic time, while their morphology showed secular changes^. In the 'admixture' model, Jomon 
admixed with the Yayoi, more recent continental immigrants, which is consistent with the rapid changes in 
morphology and culture which took place synchronically about 2,500 years before present (BP) '•'. In the 'replace- 
ment' model. Paleolithic Jomon was completely replaced by the continental immigrants (Yayoi) after their 
arrival*. To date, the 'admixture' model is seemingly better supported by the increasing lines of evidence of 
multiple genetic components found in modern Japanese^ 

The upper Paleolithic populations, i.e. Jomon, reached Japan 30,000 years ago from somewhere in Asia when 
the present Japanese Islands were connected to the continent". The separation of Japanese archipelago from the 
continent led to a long period (—13,000 - 2,300 years B.P) of isolation and independent evolution of Jomon'^. The 
patterns of intraregional craniofacial diversity in Japan suggest little effect on the genetic structure of the Jomon 
from long-term gene flow stemming from an outside source during the isolation'". The isolation was ended by 
large-scale influxes of immigrants, known as Yayoi, carrying rice farming technology and metal tools via the 
Korean Peninsula. The immigration began around 2,300 years B.P. and continued for the subsequent 1 ,000 years^. 
Based on linguistic studies, it is suggested that the immigrants were likely from Northern China, but not a branch 
of proto-Korean". 

Genetic studies on Y-chromosome and mitochondrial haplogroups disclosed more details about origins of 
modern Japanese. In Japanese, about 51.8% of paternal lineages belong to haplogroup O"", and mostly the 
subgroups 03 and 02b, both of which were frequently observed in mainland populations of East Asia, such as 
Han Chinese and Korean. Another Y haplogroup, D2, making up 35% of the Japanese male lineages, could only be 
found in Japan' '^. The haplogroups Dl, D3, and D*, the closest relatives of D2, are scattered around very specific 
regions of Asia, such as the Andaman Islands, Indonesia, Southwest China, and Tibet''\ In addition, CI is the 
other haplogroup unique to Japan"" '^. It was therefore speculated that haplogroups D2 and O may represent 
Jomon and Yayoi migrants, respectively*. However, no mitochondrial haplotypes, except M7a, that shows 
significant difference in distribution between modern Japanese and mainlanders^. Interestingly, a recent study 
of genome-wide SNPs showed that 7,003 Japanese individuals could be assigned to two differentiated clusters. 
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Hondo and Ryukyu, further supporting the notion that modern 
Japanese may be descendent of the admixture of two different 
components^. 

However, the estimation of contribution of the PaleoUthic contin- 
gency, i.e. Jomon, to modern Japanese posed a technical challenge 
since all contemporary populations in Japan were subject to admix- 
ture at various levels and the 'pure' Jomon no longer exist. The 
STRUCTURE and similar analysis does provide a solution of estim- 
ating the relative contribution of different components in admixed 
populations when both parental populations are avaUable"''^, 
but such approach was challenged for its accuracy especially when 
information on parental populations is incomplete"" ". 

In this report, we estimated the relative contribution of Paleolithic 
inhabitants (Jomon) and more recent immigrants (Yayoi) to modern 
Japanese, under the assumption of the 'admixture' model. This was 
achieved with only one of the parental populations (Yayoi) was 
available while another parental population (Jomon) was missing. 
We also explored the possibility of inferring the genetic origin of 
Jomon. 

Results 

Genetic contribution of Jomon in mordent Japanese. With the 
given two-population admixture model (Figure 1, see Materials 
and Methods section for details), the contributions of Jomon to 
Hondo Japanese (JP_ML) and Ryukyuans (JP_RK) were estimated 
by using genome-wide SNP data. Genotypes of both JP_ML and 
JP_RK populations were produced in PanAsia SNP project (see SI 
of Supplementary fUe). 

The original genetic donor to Yayoi migration is still controversial, 
both Korean and Chinese (northern or southern Han) were potential 
candidates^ '". Since genetic difference exists between northern and 
southern Han Chinese populations", we examined the possible con- 
tribution of both northern Han Chinese (NH, combining data of 
CHB in HapMap project and NHan in HGDP) and southern Han 
Chinese (SH, combining data of CN-SH in PanAsia SNP project and 
SHan in HGDP) respectively in this study. We also examined the 
contribution of northern East Asian population by combining data 
of NH and Korean (NHK, including NH and Korean population 
named KR-KR in PanAsia SNP project) since previous genetic stud- 
ies differed on whether Koreans were the only immigrants in period 
of later Jomon and early Yayoi"". Therefore totally three potential 
genetic sources (NH, SH and NHK) were evaluated respectively as 
continental genetic donor. 

The estimated contribution of Jomon to modern Hondo Japanese 
ranges from 0.231 to 0.395 for different donors (Figure 2). The 




Figure 1 Demographic model of population admixture with missing 
population. 




NHK SH NH NHK SH NH 

Figure 2 | Estimated contributions of Jomon in different genetic 
scenarios. Proportions of Jomon's contribution were marked on y-axis 
and names of the Asian mainland donors were presented on x-axis. Results 
based on parameters presented in table 1 were shown in blue columns. 
Open columns presented estimations from scenarios where Jomon's Ne is 
only a half of before. Black bars indicated 95% C.I of the estimations. 

contribution of Jomon (0.231) is the lowest, when NHK (combina- 
tion of Northern Han and Korean) was taken as the donor, with 95% 
confidential interval (C.I.) 0.215-0.266. When NH and SH were 
taken as donors, estimated contributions were 0.395 (C.l. 0.379- 
0.421) and 0.376 (C.I. 0.344-0.387), respectively. In contrast, 
Jomon's contribution in modern Ryukyuans is much greater than 
that in modern Hondo Japanese, regardless the choice of donors. The 
estimated contribution are 0.543 (95% C.I. 0.512-0.567), 0.605 (95% 
C.I. 0.605-0.648) and 0.623 (95% C.I. 0.623-0.671), for NHK, NH, 
and SH, respectively (Figure 2). 

Combining lines of evidence from the above results and previous 
reports on both maternal and paternal lineages, none of the three 
possible donors can be excluded from the model for peopling of 
Japan. We therefore averaged the above estimations from all three 
possible donors to assess the proportion of Jomon contribution in 
modern Japanese. Thus, in Hondo Japanese, 33.4% of genetic com- 
ponent was derived from Jomon, whereas, 60.5% of its genetic com- 
ponent in Ryukyuans may come from Jomon. These estimations 
based on autosome data are consistent with previous studies using 
frequencies of Y haplogoups in modern Japanese populations. In 
particular, overall frequency of Japanese specific Y lineages, D2 
and CI, is about 60% in Okinawa (Ryukyu) and 26.4—46.2% in 
Japanese mainland"". 

It is likely that the effective population size (Ne) of Jomon was 
relatively smaller than that of the populations in Asian main- 
land, because history of agriculture is relatively shorter in Japan 
Archipelago than that in Asian mainland. We applied our estimation 
in scenarios that Ne of ancestry Japanese is only a half of the size in 
aforementioned estimations. The extent evaluation shows the 
reduced Ne led to only minor changes to the estimations of Jomon 
contribution (Figure 2). In particular, the estimated contribution 
with reduced Ne is 0.194 (95% C.I. 0.146-0.220), 0.326 (95% 
C.I. 0.321-0.356) and 0.369 (95% C.I. 0.356-0.387) for Hondo 
Japanese and 0.510 (95% C.I. 0.490-0.525), 0.557 (95% C.I. 0.546- 
0.592) and 0.575 (95% C.I. 0.554-0.602) for Ryukyuans with three 
possible donors, respectively. In other words, the estimations are 
robust to the changes of the parameters of demographic model. 

Genetic affinity of Jomons. With the genetic contribution of 
Jomon in Ryukyuans and Hondo Japanese being 33.4% and 
60.5% respectively, allele frequencies of SNPs of ancestral Jomon 
populations were inferred using the maximum likelihood (ML) 
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method described in the Materials and Methods section. Overall 
allele frequencies in 6 virtual Jomon populations (JOMON ML- 
NH, JOMON ML-SH, JOMON ML-NHK, JOMON RK-NH, 
JOMON RK-SH and JOMON RK-NHK) were inferred based on 
two extant Japanese populations, i.e. Hondo and Ryukyuans, and 
three possible continental donors including NH, SH, and NHK. 
The inferred allele frequencies of Jomon populations allowed the 
reconstruction of the phylogeny including Jomon and the extant 
populations, using coancestry coefficient^" as the measurement of 
genetic distance and Neighbor- Joining method (NJ) for phylogeny 
reconstruction. Interestingly, all six inferred Jomon populations fell 
into the group of Northeast Asian populations but not that of the 
populations in South Asia or Southeast Asia (Figure 3A). The results 
of Principle Component Analysis (PCA) confirmed the genetic 
affinity shown in the NJ tree (Figure 3B). Jomon is therefore 
genetically closer to North Asian populations than it is to any 
other populations, suggesting that it is more likely of North Asian 
origin. The results based on autosomal data do not support the 
hypothesis of Southeast Asian origin which was proposed by 
Tuner 2""^ (1976)''. In addition, PCA revealed that two modern 
Japanese are located between mainland Asian populations and 
Jomon populations, consistent with the notion that they are the 
descendents of the admixture of mainland Asians and Jomons 
(Figure 3B). Furthermore, the clustering of Jomon populations 
inferred from Hondo and Ryukyuan Japanese suggested that their 
respective Jomon components were of similar, if not identical, 
origins (Figure 3A &B). 

Discussion 

Through the estimations and inferences presented above, this study 
may shed light on the understanding of peopling of Japan. A few 
critical issues related to the findings of this work need to be further 
scrutinized. 

Our result suggested D2 lineage of Y chromosome provides insuf- 
ficient information on Jomon's origin. Sharing of haplogroup D- 
M174 between central Asian populations and Japanese has been 
recognized as evidence of central Asian origin for paternal lineage 
of Paleolithic Japanese*. However, this study showed that some 
populations carrying branches of Y haplogroup D-M174 (D*: 10% 
in Dai; Dl: 7.7% in Chuang, 2.59% in Tujia; D3: 28.89% in Naxi) 
showed little genetic relationship with the inferred Jomon popula- 
tions (Figure 3)". Furthermore, Y haplogroup D-M174 appears 
rarely in North Asian populations when the populations have a close 
genetic relationship with Jomon (Figure 3). The only appearance of 
D2 lineage in Japanese and the widespread geographic distribution of 
D-M174 made it practically impossible to pinpoint a population in 
the continent that provided the D2 lineage to Japanese. In contrast, 
autosomal markers indeed provide additional information that con- 
nects Jomon and North Asians. 

In this study, a specified model was used to evaluate genetic 
admixture in modern Japanese populations. However, there is 
another powerful approach that was commonly applied in data ana- 
lysis for population admixture, Bayesian method implemented in 
STRUCTURE. The Bayesian method focuses on a genetic inher- 
itance model specified in terms of the proportion of an individual's 
genome originating from each of a set of possible subpopulations^'. 
The method is powerful and well recognized. However, when popu- 
lation history is partially understood and the divergence among 
ancestral populations is relatively small, performance of the Bay- 
esian method is not fully clarified yet^^. We applied the Bayesian 
method on 4 Japanese populations and 9 other East Asian popula- 
tions with assumption that all the populations shared 2, 3, or 4 
independent ancestry populations (Figure 4). None of inferred gen- 
etic components from STUCTURE could well explain the great 
frequency difference of Japanese specific Y lineages, D2 and CI, 
between modern Hondo Japanese and Ryukyuans, while our 



estimations were concordant with the previous reports (Figure 2 & 
3). The comparison suggested a specified model could be superior 
when population history was partially known. 

Geographic distribution of lineages explained the great contri- 
bution of Yayoi in our results. Hammer et al. investigated geographic 
distribution of Y lineages in Japanese populations. Haplogroup fre- 
quencies of the Y lineages showed U-shape cline with significant 
correlation with geographic distance of the populations from 
Kyushu. In briefs, the frequency of D2 lineage increased with increase 
of the distance meanwhile frequencies of O lineages decreased"". The 
O lineages were recognized as a Yayoi founding lineage and D2 
lineage was believed to be Jomon specific'''^''. Therefore, the pattern 
of geographic distribution of lineages supported published archeo- 
logical and anthropological results about population expansion dur- 
ing Jomon and Yayoi period in Japan. The archeological studies 
suggested general demographic density was significantly greater in 
eastern Japan compared to western Japan around the 3,300 years BP 
and a rapid increase first happened in West Japan around 2,000 years 
BP^"*. The studies of physical anthropology on human skeleton 
showed the new continental immigrants in West Japan, Yayoi peo- 
ple, have better capability to achieve enough foods to feed more 
people than Jomon^''"'^. The pattern of population expansion may 
explain the great genetic contribution (about 60-72%) of Yayoi in 
extent Japanese. Size of continental immigration was not necessary to 
be very large but descendants of the immigrants increased rapidly 
and subsequently dispersed from West Japan to other regions. Popu- 
lation admixture between the continental descendants (Yayoi) and 
Jomon descendants shaped genetic pattern of extent Japanese. Straits 
between Japanese islands and Asian mainland may not act as effec- 
tive barriers to the genetic admixture. 

This study supplied a genetic view for peopling of Japan. Both 
contributions and genetic affinity of Paleolithic Jomon were investi- 
gated in quantitative approach. However, the results should be 
improved continuously in further analysis. At first, the demographic 
model shall be improved when more data become available across 
Eurasia, especially from East Asia. Any improvement for the demo- 
graphic model will lead to better estimation and inference in the 
model-based approach. The second, the used divergence time be- 
tween historical Jomon and continental Yayoi should be validated 
with increased lines of evidence in future although both archeological 
evidences and genetic studies supported a divergence time between 
Jomon and continental populations is properly around 20,000 years 
BP^'*". Furthermore, relatives of Yayoi in Asian mainland should be 
specified further. 

Methods 

Genotype data. This study involves 94 populations from worldwide, mostly Asia. 
Genotype data of these populations were collected from three sources, Pan- Asia SNP 
project^*", the International HapMap project^^ and Human Genome Diversity Project 
(HGDP)^'. In the Pan-Asia SNP project, DNA samples from 1,719 unrelated 
individuals of 71 representative populations from China, India, Indonesia, Japan, 
Malaysia, the Philippines, Singapore, South Korea, and Thailand etc. were collected 
and genotyped in Affymetrix platform^*". Genotypes for 60 unrelated European- 
Americans (CEU), 60 unrelated Yoruba (YRI), 45 unrelated Chinese (CHB), and 44 
unrelated Japanese (JPT) were downloaded from the website of International 
HapJVIap project. Furthermore, the HGDP data was obtained from website of 
Stanford Human Genome Center for genotypes of 240 individuals of 17 East Asia 
populations and 45 individuals of 2 Native American populations. Finally, 20,362 
SNPs with average spacing of 134.8 kb shared among all the 94 populations were 
included for further analyses. Names of populations, population IDs, sample sizes and 
more details could be found in online material {see SI of Supplementary file). 

Admixture model. A genetic model (Figure 1) including two parental populations 
(Pi and P2) and one admixed population (Ph) is frequently used to evaluate 
contribution of historical genetic donors (population Pi' and P2'} to the population 
with admixture (Ph)'^^'*'^- The Pis P2' and Ph' are ancestral populations toPi,P2andPh, 
respectively at t^MD generations before present when the parental populations Pi' and 
P2' admixed and led to emergence of Ph' . Pr and P2' contributed to Ph' with the 
proportion ofX and 1 - X, respectively. AU the populations have a common ancestral 
population, Pq, at the divergence of the two parental populations, Pi' and P2'. Tune of 
the divergence is t^v generations before the admixture event. Each population in 
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Central East Asian 




Figure 3 | Genetic affinities between Jomon and other Asia populations. A. NJ tree based on pairwise genetic distance. Each inferred Jomon population 
was named in population IDs of its admixed descendants and continental donors (JOMON ML-NH, JOMON ML-SH, JOMON ML-NHK, JOMON RK- 
NH, JOMON RK-SH and JOMON RK-NHK, respectively). B. Populations presented on a 3D plot with axis of PCI, 2 and 3 from PCA analysis. Features of 
the genetic affinities were the same as those were discovered by NJ tree. 
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Figure 4 | Estimated admixtures with different pre-given value for number of genetic components in STRUCTURE 2.3. 



each historical stage have effective population sizeNj, j e {1, 2, h, 1', 2'} (see Figure 1 
for detail). Samples (in sizes of S^, S2 and Sh) were collected from the three currently 
observable populations (Pi, P2 and Ph). In Wright-Fisher model on assumption of 
Hardy-Weinberg Equilibrium (HWE), probability of the observation in samples can 
be calculated using Bayesian theorem when all demographic parameters are known.^^ 

Estimating historical genetic contribution. Likelihood method was applied on the 
aforementioned demographic model for parameter estimation^^. Wang's method, 
implemented in LEADMIX, estimates multiple demographic parameters 
simultaneously using genotype data from the parental informative populations (Pi 
and P2) and the admixed descendants (Ph). In this study, however, the genotypes on 
one of the parental populations, i.e. Jomon, is unavailable since none of the modern 
Japanese populations could be used as pure descendants of Jomon without admixture 
with Yayoi people. LEADMIX neither was evaluated for its performance for the 
scenarios with missing parental population, not it could handle thousands of markers 
simultaneously. Therefore, we developed another maximum likelihood (ML) 
approach to estimate the contribution of the missing parental population. 

The basic idea of ML method is to maximize the probability of observations with an 
appropriate value for genetic contribution X. For a specific locus j, the observed counts 
of a given allele are and Cjh in samples from population P2 and Ph, respectively. 
Given the demographic parameters {oj|52,5/i,iV2,A/;j,iVi ,/}, the probability of 
observation on locus j could be calculated as bellows 

Prfe>) f f [ [ [ ^f^2^'"2_ dfi dMn (1) 



Pr(<:|a))oc H Pr(c,|co) 



(3) 



where 



J J dfodfi' dfi dfidfh 

fo fi fj fl fh 

h(Ci\co)=PT{Ci\fi,fH,S2,Shy'Pr{f2,Mfi,fh,N2,Nh,tAMD) 

■ Prte'|;.,/i /2 )■ PrCfisMfcW, ,N2 ,tD,y)- Pr(/o) 



(2) 



Here, {f^.|k— 0, 2, h, 2', h'} is allele frequency of the given allele in a specific population 
k. Effective population sizes are given in {N[^|k— 0, 2, h, 2', h'}. 

With n independent loci or loci with ignorable linkage disequilibrium, we have 
probability of observations on the n loci 



Thus, genetic contribution X can be estimated by maximizing the probability of 
observations using the ML approach. 

During the calculation for probability of observations, allele frequencies in des- 
cendant populations were determined by allele frequency in ancestry population Po 
and other given demographic parameters, {g)|S2, Sh, N2, Nh, Ni-, N2', X}, in the 
Wright-Fisher model. The allele frequency in common ancestral population Pq could 
be determined by either coalescence theory or uniform distribution with range from 0 
to p2.29-3i However, joint frequency spectrum of given alleles in descendent popu- 
lations (Pis P2', Phs Pi> P2 ^rid Ph) is difficult to be presented in closed form. We 
utilized a forward Monte Carlo (MC) simulation to obtain the joint frequency 
spectrum and calculate the probability of observations. 

To maximize the probability of observations, an algorithm of golden section search 
was implemented with the aforementioned probability calculation. The algorithm 
finds the extremum of probability by successively narrowing the range of genetic 
contribution A.^^. In other words, we seek a genetic contribution X given the best 
explanation to the observations. We executed the search algorithm for multiple times 
on independent data sets of the same populations to reduce uncertainty of the MC 
calculation. Median of the optimized results supplies a robust estimation of genetic 
contribution X from results of the multiple runs. The two-side 95% confidential 
intervals of the median are given by order statistics^^. Evaluation on simulation data 
showed that our approach was accurate in multiple scenarios and robust to ascertain 
bias of marker allele frequency (see 82 and S3 of Supplementary file). 

Estimating contribution of Jomon. Parameters of the demographic model 
(Figure 1), except genetic contribution X, must be given before further admixture 
analysis using aforementioned ML method. 

The estimated time of divergence between Jomon and their continental relatives is 
— 14,000-30,000 years B.P., properly around 20,000 years B.P. (based on coalescent 
time of Y haplogroup M174-D2)^-^. The arrival of Yayoi people and the admixture of 
Yayoi and Jomon occurred about 2,000—2,500 years BP^. We therefore let 
^Div ~ 1,000 and t^MD ~ 100 generations in subsequent analyses. The rapid accu- 
mulation of genetic data improves the reliability of demographic model of human 
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Table 1 | Demographic parameters of the best fitted 
Time duration (generations before present) 


model for East Asia populations 

He 


Stage of He 


U--4UU 


1 r\r\c\r\c\ 
\ UUUUU 


n -11 

Post-agricultural expansion He 


400- 1 900 


7700 


Asian common He 


1900-2000 


725 


East Asian bottleneck He 


2000-3429 


7700 


Non-African He 


3429-3500 


400 


Out of Africa bottleneck Ne 


3500- 1 7000 


24000 


African He 


1 7000~before 


12500 


Human ancestral He 



populations. Based on fitness of allele frequency spectrum of simulation data and 
multi-locus sequencing data, Schaffner et al. identified an optimized demographic 
model including Han Chinese and Japanese^*. The best fitted demographic model 
with its parameters have been well supported and widely applied in genetic studies^^"'^. 
Thus, we employed the parameters of population dynamics described in Schaffner 
et al (2005) in our admixture estimation {table 1). 

In each estimation, 30 SNP sets were involved, each with 500 SNPs randomly 
selected with replacement from 20,362 SNPs. The contribution of Jomon to extant 
Japanese populations was estimated using all 30 groups individually and the median 
of the estimations was taken as the overall estimation. 

Inference of allele frequency of the missing parental population. In order to 
explore genetic affmity of Jomon with the extant populations, we estimated allele 
frequency of each SNP in Jomon population using an ML method. Given a genetic 
contribution \ and current observation Cj2 and Cjh on locus j for population P2 and 
P]^, allele frequency fji' in the missing ancestral population P^' could be estimated 
when the recent effective population sizes (Nh and N2) are large and the admixture 
event is relatively young. Details of the estimation were addressed below. 

Given the current observation Cj {Cj2, Cjh) and sample size S (82, Sh), we have the 
probability of observation 




where 

^7,'=4^-r + (1-^)^-2' 

In order to maximize the probabilityP(C/2nQ/i[^2i^/i^'52 55/,), we let partial deriva- 
tives of Equation 4 to be zero. After solving the partial derivative equations, we have 
estimation for the allele frequency in missing parental population fji by 

( Cjul\-{\-l)C^I&i Cju|\-(^-X)0,^|$^ 

h= \ ^ .^ CjH/SH-{l-3)Cj2/S2 ^ ^ (5) 

I Q^otherwise. 

Discovery of Genetic Relationships. With the estimated genetic contributions of 
Jomon in Ryukyuans and Hondo Japanese, we inferred the allele frequencies of 
corresponding markers of Jomon populations using genotype data of Ryukyuans 
(JP_RK), Hondo Japanese (JP_ML) and all the three possible genetic sources (NH, SH 
and NHK). Allele frequencies of Jomon were inferred separately utilizing pairs of the 
data (JP_RK and NH, JP_RK and SH, JP_RK and NHK, JP_ML and NH, JP_ML and 
SH, JP_ML and NHK) in this study. Using the inferred frequencies, we explored 
genetic relation between Jomon and other Asian populations. 

The coancestry coefficient is a measure of genetic distance and the coefficient 
determined by drift only is proportion to divergence time of pairwise populations^". 
We calculated the coancestry coefficient to uncover pairwise genetic distance between 
Jomon and other Asian populations. Based on the genetic distances, NJ tree was built 
to present possible genetic affinities in program Splits Tree 4.0^°. AU mordent 
Japanese populations and Korean were excluded from the NJ tree to avoid possible 
bias because coancestry coefficient was not originally proposed to handle such 
admixed populations. Simultaneously, Principle Component Analysis (PCA) was 
used to visualize genetic affinity based on direct difference of the allele frequencies of 
all populations. The PCA analysis was conducted in Matlab language and compu- 
tational environment (version 2007b). 
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